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Abstract 

The  kinetics  of  the  water  gas  shift  reaction  was  studied  to  evaluate  existing  reaction  mechanisms,  test  various  rate  expressions  and  simulate 
the  performance  in  a  methanol  fuel  processor  for  fuel  cell  applications.  The  reaction  was  carried  out  in  a  micro  reactor  testing  unit  using  a 
commercial  Sud-Chemie  Cu/ZnO/AFCb  catalyst  between  120  and  250  °C  with  a  range  of  feed  rates  and  compositions.  Using  non-linear 
least  squares  optimization,  the  parameters  in  five  rate  expressions  were  fit  to  the  experimental  data.  Based  on  a  review  of  published  work  on 
the  WGS  reaction  mechanism,  our  study  found  that  a  rate  expression  derived  from  a  regenerative  mechanism  and  another  rate  expression 
derived  from  adsorptive  mechanism  fit  the  experimental  data  equally  well.  Numerical  integration  of  a  one-dimensional  PFR  model  was 
used  for  this  parameter  fitting.  An  empirical  rate  expression,  rco  =  fePcoPmot  1  —  P)  with  activation  energy  of  47.4kJ/mol  was  also 
obtained  from  the  experimental  data.  Reactor  performance  was  simulated  to  determine  catalyst  loadings  required  to  achieve  specific  CO 
conversions  as  a  function  of  temperature  and  water  feed  rate.  These  results  are  useful  in  studying  the  design  trade  offs  available  to  reformer 
systems. 
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1.  Introduction 

The  water  gas  shift  reaction,  used  to  reduce  carbon 
monoxide  in  a  hydrogen  rich  stream,  is  an  historical  and 
industrially  important  reaction.  After  the  development  of 
the  low  temperature  shift  catalyst  composed  of  copper  and 
zinc  oxide  during  the  1960s,  the  main  application  of  the 
water  gas  shift  reaction  has  been  hydrogen  production  for 
ammonia  synthesis  or  other  industrial  processes  such  as  hy¬ 
dro  treating  of  petroleum  stocks.  However,  recently  a  new 
application  for  the  water  gas  shift  reaction  is  in  reform¬ 
ing  systems  for  fuel  cells.  In  particular,  research  activity 
is  growing  in  catalyst  characterization,  kinetic  and  reactor 
modeling  and  new  catalyst  formulation. 

To  design  an  efficient  fuel  reformer  and  to  optimize  its  op¬ 
erating  conditions,  knowing  the  kinetics  for  steam  reforming 
and  water  gas  shift  is  critical.  In  most  hydrocarbon  proces¬ 
sors,  the  water  gas  shift  reactor  is  the  biggest  and  heaviest 
component  because  the  reaction  is  relatively  slow  compared 
to  the  other  reactions  and  is  inhibited  at  higher  temperatures 
by  thermodynamics.  Therefore,  reducing  the  size  of  the  wa¬ 
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ter  gas  shift  reactor  is  an  important  issue.  For  fuel  cell  ap¬ 
plications,  a  compact,  efficient  and  reliable  fuel  processor  is 
highly  desirable.  To  design  a  compact  fuel  reformer  system 
using  process  simulation  and  optimization,  WGS  reaction 
kinetics  are  a  required  and  key  component. 

A  number  of  rate  expressions  have  been  reported  and 
tested  to  evaluate  the  water  gas  shift  reaction  rate  for  various 
catalysts  [1]: 

CO  +  H20  ++  CO2  +  H2,  AH°29S  =  -41.1kJ/mol  (1) 

Different  mechanisms,  various  elementary  reaction  paths, 
and  different  rate  determining  steps  and  assumptions  as  to 
the  active  sites  of  the  catalyst,  create  numerous  rate  expres¬ 
sions.  These  expressions  can  be  successfully  fit  to  specific 
experimental  data.  If  the  data  is  highly  precise,  the  goodness 
of  fit  can  help  guide  the  choice  of  mechanism.  For  example, 
the  rate  expression  of  a  large  industrial  WGS  reactor  op¬ 
erating  at  high  temperature  and  pressure  could  be  different 
from  a  small  reactor  in  a  methanol  reformer  for  a  PEM  fuel 
cell.  In  our  study,  we  focus  on  those  conditions  most  likely 
in  the  WGS  reaction  for  a  medium  scale  methanol  reformer 
using  a  Sud-Chemie  Cu/ZnO/AHOy  commercial  catalyst. 

Although  numerous  studies  of  the  reaction  kinetics  and 
mechanism  for  the  water  gas  shift  reaction  have  been 
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Nomenclature 

A 

constant 

Ao 

pre-exponential  factor 
(molgcat-1  h-1  atm-2) 

E 

activation  energy  (kJ/mol) 

A 

i 

standard  heat  of  reaction  (kJ/mol) 

number  of  data  points 

k 

rate  constant  for  the  water  gas  shift 
reaction  (molgcat-1  h-1  atm-2) 

k\ 

rate  constant  of  H2O  adsorption 

k2 

rate  constant  of  CO  adsorption 

k- 1,  k- 2 

rate  constant  of  the  reversible 

shift  reaction 

K,Ke  q 

equilibrium  constant 

Kp 

equilibrium  constant  of  water  gas 
shift  reaction 

Ki 

adsorption  equilibrium  constants  for  H2O 

k2 

adsorption  equilibrium  constants  for  CO 

k3 

adsorption  equilibrium  constants  for  CO2 

k4 

adsorption  equilibrium  constants  for  H2 

m 

concentration  exponent  of  H2O 

n 

concentration  exponent  of  CO 

neco.i 

calculated  molar  flow  rate  of  CO  at 
the  reactor  outlet  (mol/h) 

neCO,i 

experimental  molar  flow  rate  of  CO  at 
the  reactor  outlet  (mol/h) 

Pco 

partial  pressure  of  CO  (atm) 

Pc  o2 

partial  pressure  of  CO2  (atm) 

P\  t2 

partial  pressure  of  H2  (atm) 

/,it2o 

partial  pressure  of  H2O  (atm) 

rco 

reaction  rate  of  CO  (molgcat-1  h-1) 

R 

gas  constant  (J/mol  K) 

T 

temperature  (K) 

Greek  letter 

P 

factor  of  reversible  reaction, 

P  =  P\  h  Pah  /  Pco  Pi  ho  K P 

reported  during  the  past  few  decades,  there  are  still  dis¬ 
agreements  and  contradictory  reports  over  the  active  site 
and  the  reaction  mechanism.  It  is  also  apparent  that  there  are 
two  distinctly  different  kinetic  mechanisms:  the  “adsorptive 
mechanism”  (reaction  (2))  and  “regenerative  mechanism” 
(reactions  (3)  and  (4)): 


h2o  +  CO  - 

-*  CO(ads)  +  H20(ads)  - 

*  [int*]  ->■  C02(ads) 

T  H2(ads) 

-»  co2  +  h2 

(2) 

H2O  +  Red 

->  H2  +  Ox 

(3) 

CO  +  Ox 

CO2  +  Red 

(4) 

where  Red  represents  a  reduced  site  and  Ox  an  oxidized  site. 

In  the  adsorptive  mechanism,  CO  and  H20  adsorb  on  the 
catalyst  surface  and  form  an  intermediate  which  results  in 


desorbed  hydrogen  and  CO2.  Many  research  groups  have 
tried  to  prove  the  existence  and  form  of  the  intermediates, 
such  as  the  formate  species  through  chemical  trapping  ex¬ 
periments,  isotopic  labeling,  or  IR  spectroscopy  [2-5].  The 
regenerative  mechanism  or  redox  mechanism  is  the  cycling 
of  two  steps  (reactions  (3)  and  (4)).  In  the  first  step  water 
adsorbs  and  dissociates  on  reduced  sites  of  catalyst  surface 
to  produce  hydrogen  while  oxidizing  a  site.  In  the  following 
step  CO  is  oxidized  to  CO2  on  this  oxidized  sites  [6-8]. 

From  these  two  mechanisms,  a  variety  of  rate  expres¬ 
sions  can  be  derived.  From  the  adsorptive  mechanism, 
Langmuir-Hinshelwood  type  rate  expressions  can  be  de¬ 
rived.  The  following  rate  expression  is  derived  from 
Yang-Hougen  table  [9]  when  the  surface  reaction  is  as¬ 
sumed  rate  controlling: 

r  _  ^ _ C[i2o  Pco  -  Pn2  pco2  /  _ 

(1  +  K\Pqo  +  Ki  P]\2o  +  KjPco2  +  K4PU2)2 

(5) 


The  above  rate  expression  has  been  tested  using  plant  and 
laboratory  data  by  several  authors  [10,11]  who  report  that 
only  the  Langmuir-Hinshelwood  model  can  accommodate 
all  of  the  experimental  data.  Alternatively,  single  path  reac¬ 
tion  mechanism  applied  to  the  adsorptive  mechanism  results 
in  the  following  expressions  [12]: 


.  FhioFco  -  Pn^Pcch/Kp 

rco  —  k - 

1  +  Al  Fh,0  +  Kt,  PCo2 


(6) 


For  the  redox  mechanism,  Shchibrya  et  al.  derived  the 
following  rate  expression  and  confirmed  its  validity  using 
Cu-Zn-Cr  catalysts  [6]: 


rco  = 


OboAx)  -  Pr2  I\xhJ  K\> 
AP\uo  +  Pc  o 


(7) 


Other  research  groups  also  report  the  validity  of  this  rate 
expression  from  the  oxidation-reduction  mechanism  [13]. 
Using  a  CuO/ZnO  catalyst,  another  rate  expression  can  be 
derived  from  the  redox  mechanism  when  a  single  path  reac¬ 
tion  model  is  assumed  [12]: 

r  _  k\k2( Pyt2 o Pco  ~  P\h  PcxhJ K\>) 

k\ Fh2o  +  k2Pc  o  +  (£1-  +  k2-)Pco2 


In  contrast  to  those  rate  expressions  from  detailed  reac¬ 
tions,  mechanisms  and  rate  determining  steps,  there  are  sim¬ 
ple  empirical  rate  expressions  which  do  not  consider  any 
mechanism.  Moe  [14]  used  a  simple  reversible  rate  expres¬ 
sion  for  carbon  monoxide  conversion: 

rco  =  k  {^  PcoPwiO - C<^  H2^  =  Wco^od-fl  (9) 


where 

=  Pc  o2  Ph2 

Pco  P\  hO  Acq 

Other  research  groups  have  reported  that  the  water  gas  shift 
reaction  is  not  a  simple  order  reaction,  especially  at  higher 
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steam/CO  ratios,  and  have  tried  to  find  proper  exponent  pa¬ 
rameters  in  a  power-law  type  equation  [15]: 

rco  =  kTcoP&oQ  -  ft  CIO) 

The  water  gas  shift  reaction  is  an  important  reaction  in  the 
fuel  reformer  because  of  the  reactor’s  size.  It  is  important 
to  include  the  effects  of  other  reactions  in  the  reforming 
system:  hydrocarbon  reforming  and  CO  selective  oxidation, 
when  designing  the  WGS  reactor.  The  goal  of  this  study  is 
to  obtain  highly  accurate  kinetics  expressions  for  the  WGS 
reaction  to  create  a  tool  for  an  integrated  and  optimized 
simulation  of  a  whole  fuel  processing  system. 


2.  Experimental 

For  our  experimental  work  we  have  chosen  to  use  a 
commercial  catalyst  and  reaction  conditions  similar  to  the 
reactor  conditions  of  a  fuel  processor  such  as:  feed  gas  com¬ 
position,  space  velocity,  pressure,  and  reactor  temperature. 
The  commercial  catalyst  used  was  Cu/Zn0/Al203  catalyst 
manufactured  by  Sud-Chemie  (Catalyst  no.:  EX-2248). 
This  catalyst  is  also  good  for  methanol  decomposition  and 
reforming.  Catalyst  information  and  XPS  analysis  can  be 
found  in  our  previous  paper  [16].  The  catalyst  was  ground 
and  sieved  to  a  particle  diameter  of  200-250  (Jim  to  eliminate 
internal  diffusion  resistance  but  allow  good  gas  distribution. 

All  reaction  tests  were  performed  in  a  standard  catalyst 
test  unit.  A  stainless  steel  tubular  reactor,  1/2  in.  in  diameter 
and  12  in.  long  was  used  for  all  reaction  tests.  To  ensure 
isothermal  conditions  along  the  bed  length,  a  split  tubular 
furnace  was  used  and  the  temperature  of  catalyst  bed  was 
measured  directly  by  a  1/16  in.  J-type  thermocouple.  The 
reaction  tests  were  performed  at  temperatures  between  120 
and  250  °C.  The  water  feed  rates  were  controlled  precisely 
by  a  syringe  pump,  74900  Series  (Cole  Palmer),  from  0.5 
to  8  ml/h.  To  simulate  the  conditions  exiting  a  methanol 
reforming  unit,  the  feed  gas  stream  was  a  1:2  mixture  of 
CO  and  FF.  The  catalyst  loading  was  1.0  g  and  the  GHSV 
at  reaction  temperature  was  controlled  between  1000  and 
20,000  h”1. 

The  effluent  of  the  reactor  was  maintained  at  120  °C  with 
heating  tapes  to  avoid  liquid  condensation  and  connected  di¬ 
rectly  to  a  CARLE  Series  S  gas  chromatograph,  which  uses 
a  Hydrogen  Transfer  System  (Pd  membrane)  for  hydrogen 
analysis.  This  is  a  specially  designed  GC  has  two  thermal 
conductivity  detectors  using  two  different  carrier  gases:  he¬ 
lium  and  nitrogen.  Helium  is  the  preferred  carrier  gas  for  all 
components  except  hydrogen  while  nitrogen  is  the  proper 
carrier  gas  for  hydrogen.  Two  columns:  Alltech  Chemisorb 
107  (80-100  mesh,  6  ft  x  1/8  in.)  and  Supelco  Carboxen 
1000  (60-80  mesh,  15  ft  x  1  /8  in.)  were  connected  in  series 
to  analyze  the  condensible  and  light  gas  components.  Four 
components:  water,  H2,  CO,  CO2  were  measured  during 
each  test  run.  Material  balances  on  carbon  were  calculated 


to  verify  measurement  accuracy,  and  for  all  runs  reported 
here  were  within  3%  of  closure. 


3.  Results  and  discussion 

3.1.  Water  gas  shift  reaction  thermodynamics 

From  thermodynamic  properties  and  relations  the  equilib¬ 
rium  constant  for  the  water  gas  shift  reaction  can  be  derived 
in  a  conventional  way  as  shown  in  the  following  equation 
[17]: 

5693.5  4 

ln(Xeq)  =  +  1.077  In  T  +  5.44  x  10“4r 

-  1.125  x  10“7r2  -  -  13.148  (11) 

where 

^  Pco2pn2 
eq  Oho  !\o 

a  simpler  equation  for  Keq  is  given  by  Moe  [14]: 

/  4577.8  \ 

Keq  —  exp  (  — — - 4.33  J  (12) 

According  to  the  above  equations,  the  equilibrium  constant 
of  water  gas  shift  reaction  is  210  at  200  °C  and  38.8  at 
300  °C.  The  water  gas  shift  reaction  is  moderately  exother¬ 
mic  reaction  (Aff^g  =  — 41.1kJ/mol)  and  its  equilibrium 
constant  decreases  with  increasing  temperature. 

At  a  given  temperature  and  thus  a  given  equilibrium  con¬ 
stant,  and  a  given  feed  composition,  equilibrium  conversion 
of  the  WGS  reaction  can  be  calculated.  Fig.  1  shows  the  cal¬ 
culated  exiting  CO  concentration  if  a  reactor  reaches  equi¬ 
librium  conversion.  According  to  this  figure,  reaching  very 
low  CO  exiting  concentrations  requires  very  low  tempera¬ 
tures  and/or  high  water  to  CO  ratios. 

3.2.  Water  gas  shift  reaction  data 

Fig.  2  shows  experimental  measurements  of  carbon 
monoxide  conversion  versus  H2O/CO  ratio  over  the  tem¬ 
perature  range  of  120-250  °C.  As  expected  the  conversion 
of  CO  increases  with  increasing  H2O/CO  ratio  at  constant 
temperature.  For  example,  a  1 : 1  molar  feed  ratio  and  220  °C, 
the  conversion  reaches  70%.  The  equilibrium  conversion  for 
these  conditions  is  calculated  as  87%.  For  approximately 
an  equal-molar  H2O  and  CO  mixture,  the  effect  of  space 
velocity  was  investigated  (Fig.  3).  The  three  isotherms  in 
this  figure  show  that  as  expected  the  conversion  of  the  WGS 
reaction  decreases  slowly  with  increasing  space  velocity. 

To  determine  the  deactivation  characteristics  of  the 
Cu/ZnO/AL03  catalyst  for  the  WGS  reaction,  a  durability 
test  was  performed  for  250  h  at  various  feed  and  operating 
conditions.  For  the  case  of  methanol  decomposition  with 
this  catalyst  the  activity  decreased  rapidly  during  the  initial 
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Fig.  1.  CO  exiting  concentration  for  equilibrium  conversion  of  water  gas  shift  reaction  with  various  inlet  H2O/CO  mole  ratios,  1  atm,  feed  H2/CO  =  2.0. 


period  of  operation  [16].  However,  as  shown  in  Fig.  4,  there 
was  no  rapid  initial  decrease  of  the  activity  for  the  WGS  re¬ 
action.  The  extent  of  activity  loss  for  the  WGS  case  was  ap¬ 
proximately  10%  over  250  h  while  the  activity  of  methanol 
decomposition  dropped  almost  25%  during  the  first  10  h. 

3.3.  Kinetic  studies  cmd  rate  expressions 

To  find  the  most  accurate  rate  expression  for  the  water 
gas  shift  reaction,  several  kinetic  models  were  selected  to 
fit  the  experimental  data  of  the  30  data  sets  shown  in  the 
Fig.  2.  The  rate  expressions  selected  for  the  evaluation  are 
listed  in  Table  1 .  The  exiting  molar  flow  rate  of  each  reac¬ 
tant  and  product  was  determined  by  numerically  integrating 


a  one -dimensional  isothermal  PFR  model.  The  constants  in 
each  model  were  assumed  to  be  Arrhenius  functions  of  tem¬ 
perature  (Ao  exp(—  E/RT)). 

To  obtain  the  exiting  concentration  of  the  reactor,  the 
MATLAB  subroutine  function  ODE23  was  used  for  numer¬ 
ical  integration.  To  determine  the  kinetic  parameters,  the  ex¬ 
pression  in  Eq.  (13)  was  minimized: 

N  1 

y>eco,  -  neco,;)2  |  (13) 

(=i  J 

where  i  is  the  experiment  number  (N  =  30  for  WGS);  neco,; 
the  experimental  molar  flow  rate  of  CO  at  the  reactor  outlet 
(mol/h);  nc'co  ,  the  calculated  molar  flow  rate  of  CO  at  the 
reactor  outlet  (mol/h). 


H2O/CO  ratio 

Fig.  2.  Water  gas  shift  reaction:  CO  conversion  vs.  H2O/CO  ratio  (reaction  temperature:  175-250  °C;  catalyst  loading:  1.0  g;  pressure:  1  atm;  GHSV: 
6100  h_l). 
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Space  Velocity  (GHSV,  1/Hr) 


Fig.  3.  Water  gas  shift  reaction:  CO  conversion  as  a  function  of  space  velocity  (reaction  temperature:  155-220  °C;  catalyst  loading:  1.0  g;  pressure:  1  atm; 
H20/C0:  1.18-1.33;  GHSV:  3000-18,000  h"1). 


0  50  100  150  200  250 

Time  on  stream  (hr) 

Fig.  4.  Deactivation  of  catalyst  for  water  gas  shift  reaction  and  methanol  decomposition:  (a)  CO  conversion  vs.  time  on-stream  (reaction  temperature: 
155-220°C;  H20/CO:  1.18-1.33;  GHSV:  6100h  (b)  methanol  conversion  vs.  time  on-stream  (reaction  temperature:  300°C,  GHSV:  I  1 00 li  1 ). 


To  find  the  parameters  that  minimized  Eq.  (13),  the 
MATLAB  subroutine  function  LSQNONLIN  was  used  for 
non-linear  least  squares  optimization. 

The  quality  of  the  fit  is  demonstrated  in  Fig.  5  by  com¬ 
paring  the  observed  and  calculated  CO  exiting  flow  rates  for 


all  the  experiments  for  models  4  and  5  of  Table  1.  Among 
the  five  rate  expressions  evaluated,  model  4,  the  single  path 
redox  mechanism,  fits  the  data  best.  Model  1,  the  double 
site  Lanmuir-Hinshelwood  type  rate  expression  from  the  ad¬ 
sorptive  mechanism  also  fits  the  experimental  data  very  well. 


Table  1 

Evaluation  of  kinetic  models  from  different  mechanisms 


Model 

Type 

Rate  expression,  r 

NOPa 

Dsqb 

R 

Reference 

1 

Adsorptive 

kPlP2(l-p)/a  +  Y.liKIPi)2 

10 

0.49 

0.996 

[9,10] 

2 

Adsorptive 

kPiP2(l-  P)/(1  +  KtPt+  K3P3) 

6 

1.51 

0.989 

[12] 

3 

Regenerative 

kPlP2(\-P)/(Pi+AP2) 

4 

3.80 

0.987 

[6,13] 

4 

Regenerative 

kik2PiP2 (1  -  p)/(k1P1  +k2P2  +  K3Pi) 

6 

0.56 

0.996 

[14] 

5 

Empirical 

kPiP2(}  —  P) 

2 

1.79 

0.986 

[15] 

P\-P$:  partial  pressure  of  CO,  H2O,  CO2,  H2. 
a  Number  of  parameters. 
b°sq  =  Etifne,  -nc,)2. 
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Fig.  5.  Comparison  of  experimental  exiting  CO  flow  rate  data  with  cal¬ 
culated  numerical  value  using  models  4  and  5  in  Table  1  (temperature: 
120-250  °C;  catalyst  loading:  1.0  g;  pressure:  1  atm). 


However,  all  five  mechanisms  gave  very  high  R2  values  in¬ 
dicating  that  any  of  the  models  are  adequate  for  simulation 
purposes. 

The  empirical  rate  expression  derived  from  the  numerical 
fitting  is: 


rco  —  2.96 x  105  exp 


47400  \ 
RT  ) 


Pco2  Ah  \ 
Ke  ) 


(14) 


The  activation  energy  for  the  empirical  model  of  47.4  kJ/mol 
is  consistent  with  other  values  in  the  literature.  The  compar¬ 
ison  of  activation  energy  and  frequency  factors  with  other 
first  order  rate  expressions  or  empirical  power-law  equations 
are  summarized  in  Table  2. 


3.4.  Simulation  of  the  water  gas  shift  reactor 

The  design  of  a  water  gas  shift  reactor  used  to  generate 
fuel  cell  grade  hydrogen  must  choose  several  operating  pa¬ 
rameters  subject  to  several  constraints.  Five  critical  design 
constraints  include  (1)  the  inlet  gas  composition;  (2)  the  in- 


Table  2 

Parameter  comparison  for  empirical  expressions  (r  =  ko  exp(— E / RT) 

PcoPw od  -  W 


Catalyst 

m 

n 

ln*0 

E  (kJ/mol) 

Reference 

ICI-Cu/Zn0/Al203 

1 

1 

15.2 

52.8 

[15] 

Cu/A1203 

1 

1.9 

- 

69.3 

[8] 

Cu/Zn0/Al203 

0 

1 

- 

41.8 

[18] 

CuO/Mn02 

1 

1 

- 

55.0 

[19] 

Cu/ZnO/Al203 

1 

1 

12.6 

47.4 

This  work 

Table  3 

Design  basis  for  WGS  reactor  simulation  and  1  kW  fuel  cell 


Design  variables 

Values 

WGS  reactor 

CO  inlet  molar  flow  rate  (gmol/h) 

11.1 

Fb  inlet  molar  flow  rate  (gmol/h) 

22.2 

CO  target  conditions  (mol%) 

1 

Fuel  processor 

Required  amount  of  hydrogen  (gmol/h) 

33.3 

12.4“ 

Fuel  cell  stack 

Target  power  (kW) 

1 

Stack  voltage,  50  cells  (V) 

35 

Cell  amperage  (A) 

28.57 

Fuel  cell  efficiency,  assumed  (%) 

80 

a  The  unit  of  the  value  is  slpm. 


let  gas  flow  rate;  (3)  the  temperature  exiting  the  reforming 
reactor;  (4)  the  exiting  CO  concentration  constraint  of  the 
down  stream  processes;  and  (5)  the  water  requirements  for 
proper  operation  of  the  fuel  cell.  For  the  case  of  methanol 
steam  reforming,  the  gas  typically  contains  H2  and  CO  in  a 
ratio  between  2  and  3,  with  some  residual  water.  The  exit¬ 
ing  reformer  temperature  is  near  300  °C,  and  the  flow  rate 
must  be  high  enough  to  meet  the  fuel  cell  electrical  output 
requirements,  which  is  approximately  121  of  H2/(minkW) 
of  electrical  output  (Table  3).  Exiting  the  WGS  reactor,  the 
CO  concentration  must  be  low  enough  to  be  fed  directly  to 
the  fuel  cell  (about  50  ppm  of  CO),  or  more  typically  it  must 
be  moderately  low  so  it  can  be  treated  by  preferential  ox¬ 
idation  (about  1  mol%).  The  water  content  should  be  near 
12mol%  for  proper  humidification  of  a  PEM  fuel  cell.  The 
actual  optimum  water  concentration  depends  on  fuel  cell  op¬ 
erating  pressure  and  temperature.  The  values  of  these  con¬ 
straints  are  somewhat  flexible  since  the  full  optimization  of 
the  fuel  cell  design  could  include  them  as  independent,  not 
fixed  variables. 

The  WGS  reactor  design  parameters  used  to  meet  these 
constraints  include  (1)  reactor  temperature;  (2)  water  addi¬ 
tion  rate;  and  (3)  reactor  size  (weight  of  catalyst).  This  multi 
constrained,  three-dimensional  problem  lends  itself  to  sim¬ 
ulation  for  analysis,  because  the  tradeoffs  between  the  pa¬ 
rameters  are  not  obvious. 

For  simulation  and  analysis,  an  electrical  generating  ca¬ 
pacity  of  1  kW  is  chosen,  and  methanol  is  chosen  as  the 
reformer  feed.  The  basis  for  this  fuel  cell  is  summarized  in 
the  Table  3.  The  simulation  uses  the  empirical  kinetic  model 
shown  in  Eq.  (14)  because  of  its  accuracy  and  simplicity. 
It  also  assumes  no  mass  or  heat  transfer  resistance  and  it 
assumes  the  reactor  is  plug  flow  and  isothermal.  Although 
these  assumptions  are  not  entirely  correct,  they  allow  obser¬ 
vation  of  the  general  reactor  behavior. 

Fig.  6  is  a  contour  plot  of  the  variable  space  needed  (re¬ 
actor  temperature,  catalyst  weight,  and  feed  water  rate)  to 
meet  the  conditions  of  1  mol%  CO  exiting  the  WGS  re¬ 
actor.  Fig.  7  plots  results  from  the  same  set  of  simula¬ 
tions  as  in  Fig.  6,  except  the  vertical  axis  is  the  water  con¬ 
centration  leaving  the  WGS  reactor.  Due  to  kinetic  limi- 
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Fig.  6.  Simulation  result  of  water  gas  shift  reactor  for  1  kW  fuel  cell:  required  water  amount  to  achieve  1  mol%  of  exiting  CO  concentration  according 
to  various  temperature  and  catalyst  loading  (temperature:  200-300 °C;  CO  inlet:  ll.lmol/h;  pressure:  1  atm). 
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Fig.  7.  Simulation  result  of  water  gas  shift  reactor  for  1  kW  fuel  cell:  water  outlet  concentration  to  achieve  1  mol%  of  exiting  CO  concentration  according 
to  various  temperature  and  catalyst  loading  (temperature:  200-300 °C;  CO  inlet:  ll.lmol/h;  pressure:  1  atm). 


tations,  these  figures  show  a  steep  increase  in  water  re¬ 
quired  and  water  exiting  at  low  temperatures  and  low  cat¬ 
alyst  loadings.  Due  to  equilibrium  constraints,  they  show  a 
slow  increase  in  water  exiting  and  water  required  at  higher 
temperatures.  A  major  concern  for  fuel  cell  operation  is  to 
maintain  the  proper  water  concentration  in  the  feed  gas. 
This  suggests  it  may  be  optimal  from  a  control  perspec¬ 
tive,  to  operate  the  WGS  reactor  at  a  temperature  where  the 
water  concentration  variation  is  small  for  a  given  pertur¬ 
bation  in  feed  composition,  temperature  or  space  velocity. 
For  this  catalyst  that  temperature  would  be  at  the  bottom 


of  the  valley  in  Figs.  6  and  7,  or  near  a  temperature  of 
230-240  °C. 


4.  Conclusions 

Several  kinetic  rate  expressions  can  adequately  model  the 
water  gas  shift  reaction  over  the  Cu/ZnO/AhCT  catalyst  used 
in  this  study.  Two  rate  expressions,  one  derived  from  an 
adsorptive  mechanism  and  one  from  a  redox  mechanism, 
were  found  to  be  the  most  accurate  models.  However,  even 
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an  empirical  rate  expression  (Eq.  (14))  was  found  to  fit  the 
data  with  a  high  degree  of  accuracy. 

Applying  the  empirical  kinetic  relationship  in  a  simple 
isothermal  plug  flow  model  shows  that  to  maintain  a  con¬ 
stant  CO  and  water  concentration  exiting  the  reactor  requires 
controlling  the  reactor  temperature  and  the  water  addition 
rate  closely.  These  two  exiting  compositions  are  important 
to  maintain  steady  operation  of  the  fuel  cell.  Large  varia¬ 
tions  in  CO  concentration  will  make  downstream  oxidation 
of  CO  more  difficult  to  control,  resulting  in  fluctuating  CO 
concentrations  entering  the  fuel  cell,  which  will  result  in 
fluctuating  current  and  voltages.  Large  variations  in  water 
exiting  the  WGS  reactor  will  make  controlling  the  moisture 
needed  for  proper  PEM  operation  more  difficult. 

The  results  of  this  study  are  expected  to  be  an  important 
component  for  the  overall  design,  optimization  and  control 
of  commercial  fuel  reformers  used  to  generate  hydrogen  for 
fuel  cells. 
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